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ABSTRACT 

We present results from two statistical analyses applied to an neutral hydrogen (HI) dataset of the 
nearby tidal bridge in the Magellanic System. Primarily, analyses of the Spatial Power Spectrum 
suggest that the Magellanic Bridge, historically considered to be a single contiguous feature, may in 
fact be a projection of two kincmatically and morphologically distinct structures. The southern and 
more obviously turbulent parts appear to show structure organized similarly to the adjacent Small 
Magellanic Cloud (SMC) , while the northern regions are shown to be relatively deficient in large scale 
power. The extent of modification to the spatial power index by the velocity fluctuations is also highly 
variant across these parts of the Bridge. We find again that the northern part appears distinct from 
the southern part and from the SMC, in that the power spectrum is significantly more affected by 
slower velocity perturbations. 

We also probe the rate of spectral variation of the Hi by measuring the Spectral Correlation Function 
over selected regions. The results from this analysis highlight a tendency for the Hi spectra within 
the bright parts of the Bridge to have a more persistent correlation in the E-W direction than in the 
N-S direction. These results are considered to be quantitative evidence for the tidal processes which 
are thought to have been active throughout the evolution of the Magellanic Bridge. 
Subject headings: ISM: structure — galaxies: interactions — Magellanic Clouds — ISM: atoms 



1. INTRODUCTION 

The Magellanic Bridge (MB) is an elongated, pre- 
dominantly neutral hydrogen (Hi) structure, located be- 
tween and joining the Large and Small Magellanic Clouds 
(Kerr, Hindman, & Robinson, 1954). Numerical simula- 
tions (Ruzicka, 2003; Gardiner, Sawa, & Fujimoto, 1994; 
Gardiner & Noguchi, 1996) and the relative geometrical 
arrangement of the components of the Magellanic Sys- 
tem (e.g. Mathewson, Schwarz & Murray, 1977) indi- 
cate that the MB was formed directly from the Magel- 
lanic Clouds, following a particularly close mutual pass 
^200 Myr ago. Recent high resolution observations have 
shown that the Hi in the MB has a complex and chaotic 
filamentary structure (Muller et al. 2003). Apparently 
coherent Hi features are identifiable across the entire 
range of observed angular scales, from ~98" up to ~7°, 
although the filamentary component is clearly unevenly 
distributed throughout the entire MB region. 

Statistical techniques are a useful means to compare 
populations of similar objects in different systems, under- 
stand and model general properties and behaviors, and 
expose and isolate any hidden underlying trends within a 
dataset. In particular, statistical analyses provide some 
of the most useful descriptions of turbulent flows. The 
results of statistical analyses of emission-line data sets 
can be directly compared with theoretical predictions to 
determine the properties of the turbulence under that 
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theoretical framework. While there are many statistical 
techniques that have been used to analyze different data 
sets, this paper focuses on two complimentary methods: 
the Spatial Power Spectrum (SPS) and the Spectral Cor- 
relation Function (Rosolowsky et al, 1999). 

The SPS and the SCF were chosen because they high- 
light the importance of structure across a wide range 
of spatial scales. The organisation of spatial scales has 
a direct influence on the efficiency of energy transport 
throughout the interstellar medium (ISM) of the subject 
system and also vice versa: the transport of energy is an 
important operator on the arrangement and distribution 
of spatial scales. 

The relationship between the observed SPS and the 
characteristics of the underlying turbulent flows has 
been the work of several theoretical studies (Lazarian 
& Pogosyan 2000; Goldman 2000; Miville-Deschenes et 
al. 2003). In contrast with the SPS which measures 
structure of data sets on a more global scale, the SCF 
measures local variations as a function of spatial scale. 

No studies have yet made a direct and dedicated com- 
parison of the outcomes and predictions of the SPS and 
the SCF. Measurements of the SCF on the Hi of the LMC 
(Padoan et al. 2001) were shown to be in rough, though 
not particularly good agreement with earlier studies us- 
ing the SPS (Elmegreen et al. 2001) and a more detailed 
discussion of the predictions of the SCF in the context of 
the turbulent and small-scale structure of the LMC Hi 
were not entered into. 
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To date, studies of the SPS have been made of Hi in 
the SMC and the LMC (Stanimirovic et al. 1999a; Stan- 
imirovic & Lazarian 2001; Elmegreen et al. 2001). SCF 
studies have been made only of the LMC (Padoan et al. 
2001). These results have highlighted the existence of 
important active processes and structure trends within 
these systems and have been used as a quantifier for the 
organisation of structure within these systems indepen- 
dently. The MB forms a physical link between the SMC 
and the LMC, and represents a feature which has evolved 
purely in response to tidal forces. The MB is therefore 
an ideal system to test and compare the response of the 
SCF to a) a system which has developed in response to 
large-scale forces and b) in the context of the host sys- 
tems, i.e. the SMC and LMC. Such an analysis may 
characterize the turbulence in tidal features in general 
and in addition, a comparison of the SPS results derived 
for different components in the Magellanic System can 
provide insights into different evolutionary paths for or- 
ganization of Hi structures in an galactic versus an in- 
tergalactic medium. 

The structure of this paper is organized as follows. 
In Section 2 we review the SPS and SCF. In Section 
we summarize the observations of the subject Hi MB 
dataset. Sections 0] and [5] detail the methodology and 
results of the calculation of the Hi SPS and SCF on the 
MB Hi dataset. Sections presents the discussion and 
interpretation of our results. Finally, we summarize our 
results in Section [7| 

2. STATISTICAL TECHNIQUES 
2.1. The Spatial Power Spectrum 

The SPS has traditionally been applied to obtain a 
useful insight into hydrodynamics, MHD theory and tur- 
bulent motions (e.g. Lazarian & Pogosyan 2000). For 
example, a shallow SPS results from a system which is 
dominated by small-scale structure, while a steep SPS 
results from a system which is dominated by large-scale 
fluctuations. The SPS has been commonly applied on 
large datasets in an attempt to parameterise the struc- 
ture hierarchy of the Hi in these systems. (Crovisier & 
Dickey 1983; Green 1993; Stanimirovic et al. 1999a; Sta- 
nimirovic & Lazarian 2001; Elmegreen, Kim, & Staveley- 
Smith 2001; Dickey et al. 2001). 

The SPS is derived from the Fourier transform power 
spectrum of an observed 2-d image, denoted P(k), where 
k is the spatial frequency. Different types of observations 
commonly show that the SPS of intensity fluctuations 
can be characterized by a power-law, or the sum of two 
power laws: 

P(k) oc fc 7 or P(k) oc GW 1 + C 2 fc 72 . 

The composite power law is observed when the dynam- 
ics of a system are shaped by a process (or processes) 
which generate excess power over a specific range of 
scales. For example, supernovae are capable of inject- 
ing power across scales of tens of parsecs. The constants 
C\ and Ci are dependent on the relative importance of 
the two terms (see also Elmegreen & Scalo, 2004; Scalo 
& Elmegreen, 2004, for a detailed review of turbulent 
processes). 

In the case of a Kolmogorov-type of turbulent fluid 
(isotropic and incompressible turbulence, characterized 
by energy injection on large scales with a dissipationless 



cascade of energy over intermediate scales and dissipa- 
tion of energy only on smallest scales) the SPS of both 
3-d density and velocity fluctuations is a power law with 
a slope equal to —11/3. 

However, a spectral index of 7 ~ —3.7 does not nec- 
essarily imply Kolmogorov-like turbulence. Recent work 
using simulated datasets have shown that other turbu- 
lent regimes are capable of producing spectral indices 
which are approximately equal to this value. For ex- 
ample, Miville-Deschenes, Levrier & Falgarone (2003) 
have developed synthetic fractal Brownian motion (fBm) 
datasets which can yield spatial spectral indices of —3.7. 
Furthermore, synthetic datasets containing randomly 
distributed, supernovae-driven expanding shells also ap- 
pear to produce a Kolmogorov spectrum (Hodge & Desh- 
pande, in preparation). Both of these studies are in need 
of further clarification however, specifically in that the 
fBm set does not represent coherent turbulent structures 
(i.e. it is constructed from Brownian noise), and that the 
Hodge & Deshpande analysis involves unphysical forma- 
tion processes (i.e. they do not solve the hydrodynamic 
equations in the development of the set) . To distinguish 
and disentangle the nature and extent of velocity fluc- 
tuations from static hierarchical density structure, tech- 
niques which measure the velocity statistics of the set 
must be used. 

An important advantage the SPS has over other sta- 
tistical techniques is in relating statistics of the observed 
intensity fluctuations with the statistics of the underly- 
ing 3-d density and velocity fields. A recent development 
in SPS studies involves the averaging of data across a 
number of observed velocity channels. By integrating a 
spectral-line data cube over varying ranges of velocity, 
Lazarian & Pogosyan (2000) showed that it is possible 
to disentangle the influences of density and velocity fluc- 
tuations in the turbulent flow, i.e. the power-law index 
is a function of both the velocity, and the velocity width 
of the observations: 7 = 7(1;, Av). When an image is 
generated from a small range in velocities (small Av), 
a significant fraction of the observed structure results 
from the velocities of the turbulent flow moving emis- 
sion into different velocity channels. Since an SPS is 
usually inferred to characterize the density fluctuations 
in a turbulent flow, the velocity fluctuations would cor- 
rupt the results. However, when integrated over a large 
velocity range (large Av), most of the velocity fluctu- 
ations average out, leaving only static density fluctua- 
tions. By measuring the power law spectrum over a range 
of widths in the velocity window, Lazarian & Pogosyan 
(2000) showed that it was possible to separate the power- 
law spectrum of the density fluctuations from those of the 
velocity fluctuations. In their Velocity Channel Analy- 
sis (VCA), they refer to the regime where the velocity 
thickness of the integrated slice averages out the velocity 
fluctuations as the "thick slice regime" and the opposite 
regime where the velocity fluctuations are important as 
the "thin slice regime" . 

2.2. The Spectral Correlation Function 

As the SPS is derived from the Fourier transform of the 
data set, information from the entire image contributes 
to the power spectrum. However, this global approach 
compromises the usefulness of the SPS in isolating local 
features in the data. 
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By contrast, the SCF does not transform the data and 
operates on entire spectra in a data cube instead of en- 
tire images. The SCF parameterizes the 'similarity', So, 
between two spectra as a function of their spatial sep- 
aration, Ar and often shows a power-law relationship 
So oc Ar a over a range of separations (e.g. Padoan, 
Rosolowsky & Goodman, 2001). The SCF is normalized 
so that a perfect match between two spectra will yield a 
value of Sq = 1, and a value of zero is obtained for no 
similarity, or for an anti-correlation. Since this technique 
does not use a Fourier transform it does not require the 
data to be edge tapered and thereby reduces problems 
arising from the finite sampling and the associated edge 
effects. Therefore, the minimum size of the region that 
may be sampled with this technique is much smaller than 
that available to the power spectrum analysis. However, 
the SCF is sensitive to the noise in the data cube and its 
application is limited to high signal-to-noise data. 

The early forms of the SCF technique were used to de- 
rive an estimate of spectral similarity over a lag range 
which was defined by the velocity dispersion of the spec- 
tra of interest. It was therefore limited to data cubes 
that were described by a simple, single velocity compo- 
nent. Later modifications to the SCF made by Padoan, 
Rosolowsky & Goodman (2001) included a parameter 
which weights the function by the ratio of the RMS of 
the signal and noise. These modifications improved the 
ability of the algorithm to handle noise and expanded 
its application to more complex and multi-component 
spectra and across any range of spatial lags (Padoan, 
Rosolowsky & Goodman 2001). An interesting applica- 
tion of the SCF was shown in the study by Padoan et al. 
(2001) where a break in the power-law behavior of the 
SCF was interpreted as a signature of the LMC's disk 
thickness. 

As the SCF is a rather recent development, the re- 
lationship between the SCF results and the proper- 
ties of the observed turbulence are not understood as 
completely as for the SPS. We have attempted here a 
more dedicated study and comparison of the the SCF 
in the context of predictions from the relatively well- 
parameterized SPS. 

3. THE MB Hi DATASET 

We use the high-resolution Hi observations of Muller et 
al. (2003) for the present statistical investigation of the 
Hi distribution in the MB. These observations encom- 
pass a ~ 7° x 7° region, centered on RA = 2 h 08 m 40 s , 
Dec. = -73°27'59" (J2000), and covering the heliocen- 
tric velocity range 100-350 km s _1 . The data were col- 
lected with both the Australia Telescope Compact Ar- 
ray (ATCA) and the Parkes Telescope 1 . The observing 
parameters, data reduction and data combination pro- 
cedures are described in Muller et al. (2003). The fi- 
nal Hi data cube has l-cr brightness temperature sensi- 
tivity of ~0.9 K (corresponding to a column density of 
~1.62xl0 18 cm -2 ), a velocity resolution of ~1.6km s , 
and an angular resolution of 98". These observations are 
sensitive to all spatial scales from about 29 pc (corre- 
sponding to angular scale of ~98") to 7.3 kpc (corre- 

1 The Australia Telescope Compact Array and Parkes telescopes 
are part of the Australia Telescope which is funded by the Com- 
monwealth of Australia for operation as a National Facility man- 
aged by CSIRO. 



sponding to angular scale of ~7°), assuming a distance 
of 60 kpc to the MB. An important note in the context of 
this work, is that the ATCA interferometer, being com- 
posed of discrete elements is not capable of uniformly 
sampling the entire UV plane. The missing data is ex- 
trapolated during the deconvolution and inversion pro- 
cess which also homogonises the sensitivities somewhat 
(refer to Muller et al. 2003 and also Stanimirovic, 1999 
for explicit information on the data deconvolution and 
inversion process). 

Fig-Hshows an RA-Dec. peak intensity image of the Hi 
spectral line data cube used in this study. The diverse 
range of morphologies in the Hi of the MB is striking 
in this image: some regions are bright and filamentary, 
particularly in the southern parts, while others are more 
tenuous and smooth. Fig. shows a peak pixel image 
of the same dataset viewed in the Dec. -Velocity projec- 
tion. The distribution of Hi in this figure appears well- 
organised into at least two components. Notably, these 
two components appear to be distinct also in declina- 
tion. The velocity separation of these components is ~ 
40 km s -1 . 

The rapid spatial and velocity variability of the Hi in 
the MB shown in Figs. 13 indicates a turbulent and 
highly inhomogcncous morphology over the entire range 
of observed spatial and velocity scales. Consequently, we 
have applied a separate SPS analysis on four subregions 
within the data set, as shown in Fig. ^ with the overlaid 
grid. Each region spans ~ 2.2 x 2.2 kpc (256x256 pixels). 
We will refer to these regions as (clockwise starting from 
the bottom left on Fig. [TJ: the South East (SE), the 
North East (NE), the North West (NW) and the South 
West (SW). 

The SW region is dominated by bright, filamentary 
features across its entire extent. The Hi in the SE region 
is similarly filamentary, although the spatial coverage of 
the brightest feature is less complete. The NE region 
defines an area of less diverse structure where the Hi ap- 
pears relatively tenuous and smooth. The NW region 
appears to be dominated by part of a large loop-like fea- 
ture. Fig.|2|shows that the northern two sub-regions (NE 
and NW) encompass Hi which occurs at lower Heliocen- 
tric velocities. The SPS analyses was applied on each 
region separately. 

Due to the unpredictable response of the SCF to noisy 
data, we have restricted the SCF analysis to a the bright- 
est region of the Magellanic Bridge, shown in Fig. |3 
Much of this a rea has a spectral quality of Q > 6 (see 
Eq.[2]in Section l5~TJl . Restricting the tested data to those 
with a quality of this level limits the effects of noise on 
the SCF to < 20%. In contrast, the northern regions of 
the data set is more dominated by noise, having areas 
where Q < 6. Although binning can increase the spec- 
tral quality, the degree of averaging necessary to effect a 
substantial increase in Q ultimately obscures the velocity 
structure, thereby compromising the utility of the SCF. 
The 10x3 grid overlaid on Fig. shows the eight areas 
over which the SCF was averaged (see Each small 
square defines a 500 x500 pc area (63 x 63 pixels). As 
the telescope beam subtends approximately three image 
pixels, a new data cube was generated that sampled ev- 
ery third pixel in the x and y directions to ensure that 
the SCF samples independent spectra. 
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Fig. 1. — Peak intensity map of the entire MB Hi dataset. The overlaid square indicates the four 256 X 256 pixel (~2.2x2.2 kpc) 
subsections over which the power spectrum analysis is calculated. 



4. THE SPATIAL POWER SPECTRUM ANALYSES 
4.1. Derivation of the SPS 

The two-dimensional SPS used here is derived fol- 
lowing methodologies used in previous studies (Cro- 
visier & Dickey 1983; Green 1993; Stanimirovic et al. 
1999a; Stanimirovic & Lazarian 2001; Elmegreen, Kim, 
& Staveley-Smith 2001; Dickey et al. 2001). We start 
with two-dimensional images of the Hi intensity distribu- 
tion. These images are tapered with a Gaussian function 
to suppress the sin(fc)/fc ringing that appears from the 
sharp edges of the original image under the Fourier trans- 
form. The power (square of the modulus of the Fourier 
transform) is then measured as a function of wave num- 
ber (or spatial scale). The power is then azimuthally 
averaged in rings of logarithmically increasing intervals 
of wave numbers. 

The interferometric (ATCA) part of this dataset has 
negligible non-circularity of the u — v sampling function 
(see also Muller et al. 2003, Fig. 1) and no compensa- 
tion was made for any ellipticity of the u — v sampling 



function. The median Hi column density for the final Hi 
data cube is quite low (^4xl0 20 cm -2 ), so any correc- 
tions for self-absorption are unnecessary (see discussion 
by Minter 2002). 

We emphasise that this dataset is a combination of 
observations obtained with the Parkes single-dish tele- 
scope and with the ATCA interferometer. Therefore, a 
complete range of spatial frequencies was sampled from 
zero spacing, up to ~ 600A, corresponding to a minimum 
spatial scale of ~ 100 pc (note that this is the contigu- 
ously sampled range of spatial frequencies; although spa- 
tial frequencies up to ~ 1600A have been sampled, these 
are not part of a complete contiguous spatial frequency 
distribution). 

A sample SPS from the MB dataset, typical for velocity 
channels with significant Hi emission, is shown in Fig.^J 
Data points calculated from a continuous range of spatial 
frequencies < 600A were used for fitting and are denoted 
by squares. The error bars in the figure correspond to 1- 
a statistical uncertainties, which were taken into account 
during the fitting procedure. 
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Fig. 2. — Integrated intensity map of the Declination- Velocity projection of the combined ATCA-Parkes data cube. The northern part 
of the Hi in the Bridge is clearly offset to higher velocities with respect to the southern part by ~40-50 km s _1 . We see a clear bimodal 
distribution of Hi in the more southern part in this figure. This apparent two-fold structure originates in the body of the SMC and extends 
approximately 4 kpc towards the LMC (Muller et al. 2003). It has clearly been significant in the evolution of the Bridge and is the subject 
of a future work. The grey scale is a linear transfer function as shown in the wedge and units are K arcsec . As this map is integrated 
along the x image axis, the labeled Declination axis is correct only for RA=2 h 08 m 8 s . 
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Fig. 3. — Peak pixel map of a subregion of the MB Hi dataset. The overlaid grid shows the eight regions for which the 5(Ar) maps 
are calculated. Each square labeled with a grid reference forms the centre of a larger square region over which the SCF is calculated. The 
surrounding regions constitute the buffer. This entire area approximately corresponds to regions SE and SW of Fig. HI 



At spatial frequencies >700A the effect of the ATCA 
beam becomes prominent, as shown with a dashed line 
in Fig. The SPS derived from an emission- free velocity 
channel is shown in the same figure with crosses. We see 
that the contribution of power to the SPS from the noise 
component of the Hi dataset is negligible at low spatial 
frequencies (by a factor of 10 5 /±0 6 ' 5 «0.03 at «700A) 
and is important only at the largest spatial frequencies 
(> 700 A). As these frequencies are usually compromised 
by the ATCA beam, they are already excluded from the 
fit. Although some efforts have been made to compen- 
sate for the effects of noise and beam shape in power 
spectra (Green 1993), we have minimized these effects 
by restricting the spatial frequency range to < 700A from 



the fit. 

4.2. The SPS as a function of velocity 

Fig. [5] shows a typical SPS for all four sub-regions, 
derived from single velocity channels with significant Hi 
emission. A single power- law function, P{k) oc fc 7 , can 
be fitted to the SPS for most of the spatial frequency 
range. Similarly featureless SPS were found for the SMC 
(e.g. Stanimirovic et al. 1999a, b) and for the Galaxy 
(Green 1993; Dickey et al. 2001). 

The variation of the power-law slope with heliocentric 
velocity j(v, Av = 1.6 km s _1 ) for each of the four re- 
gions defined in Figgis shown in Fig. El The left ax is 
corresponds to the power-law index 7(1;), while the right 
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Fig. 4. — Sample SPS from the South East sub-region (dotted line), at velocity ~169 km s — 1 (heliocentric). The spectrum is typical for 
channels which contain bright Hi. The fitted line shown (solid line) is an error-weighted fit, derived using the points marked with squares, 
fmportantly, we see that the power spectra of the Bridge are able to be well fitted with a single component power-law. Also shown on this 
plot is the power spectrum for a spectral channel containing only noise (crosses) and for the 98" beam (dashed line). The error bars mark 
one standard deviation, those showing a downward pointing arrow indicate a large lower error limit. 



axis shows the mean integrated Hi brightness tempera- 
ture for a given velocity channel. 

The power-law index generally varies significantly with 
velocity, even over the velocity channels with the bright- 
est Hi emission. The nature of this variation is different 
in the different quadrants. This large variation in j(v) 
contrasts to that which was found for the Hi in the rel- 
atively homogeneous SMC: Stanimirovic et al. (1999a) 
measured no dramatic or large changes in the power-law 
index across the entire (~200 km s _1 ) velocity range. 
This has important ramifications in the next part of the 
SPS analysis, where the fundamental assumption is that 
the structure of the set under study is largely homoge- 
neous throughout the velocity range of the data. 

To some extent, the plots of Fig. suggest a relation- 
ship between the mean brightness temperature and the 
index of the power-law, in that a hight mean T^ is as- 
sociated with a steeper power spectrum. This is can be 
seen for velocity channels where the mean Tb <2K or so. 
This effect is primarily a manifestation of the increasing 
relative importance of the noise component: for channels 
where Tb is low, the spectrum becomes 'whiter' (i.e. the 
spectrum flattens). 

In the case of the SW region, the variation of T# is 
very large and we see that the variations of T^ are not 
reflected by similar variations of 7, which appears to have 
a limiting slope of approximately —3. This is close to 
that expected from a Kolmogorov flow. 



4.3. The SPS as a function of velocity slice thickness 

We now investigate the extent to which velocity fluc- 
tuations affect the apparent SPS of the Hi intensity dis- 
tribution by applying the VCA technique (Lazarian & 
Pogosyan 2000) to Hi data cube subsets averaged over 
a range of velocity intervals, which is just the measure- 
ment of j(Av). As the VCA technique assumes that the 
morphology of the original dataset is quite homogeneous 
over the velocity range of interest, we will apply the VCA 
analysis only over selected contiguous velocity windows 
for which the SPS index is approximately constant to fa- 
cilitate comparison with theoretical models. In doing so, 
we avoid mixing kinematically different regions. Unfor- 
tunately, choosing velocity windows limits the available 
velocity thickness range considerably. 

Velocity windows used for this analysis are shown as 
shaded bands in Fig. [S] We selected regions with roughly 
constant 7(1?, Av = 1.6 km s _1 ) and mean Hi brightness 
temperature (> 4 K) as a sample of relatively homoge- 
neous regions. There is no usefully large velocity range 
over which the power-law index is particularly constant 
in the NW region, so we do not attempt to measure 
7 (Aw) in this subregion. 

Quantitatively, the ranges selected for each region 
follow: 

North West: None 

North East: 196 - 230 km s" 1 

South West: 143 - 203 km s _1 
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Fig. 5. — Sample SPS from each of the four regions (marked in upper left) shown in Fig.Q These spectra are typically well fitted by 
a single component power-law. Error bars indicate one standard deviation and those showing downward pointing arrows indicate a large 
lower error limit. 



South East: 156 - 194 km s" 1 

The variation of the SPS power-law index as a function 
of velocity slice thickness for the three regions is shown in 
Fig. Generally, 7(Aw) becomes steeper with increas- 
ing velocity thickness for all three regions. This trend 
agrees with predictions by Lazarian & Pogosyan (2000) 
and suggests that both density and velocity fluctuations 
contribute to the measured Hi intensity fluctuations for 
thin slices. The selected range of velocity channels for 
the SE region is narrower than that of the SW region; 
however over the common range, j(Av) is very similar. 
This suggests that turbulence has comparable character- 
istics in the two southern regions although, the power-law 
index derived for the NE region appears very different. 
Over the whole range of Av the slope for the NE region 
is significantly shallower than the slope of the southern 
regions. 

It is interesting to note that for all three regions j(Av) 
starts to converge a constant value around Av ~ 10 
km s _1 . Lazarian & Pogosyan (2000) predicted that 
the presence of gas at different temperatures will show 
a characteristic 'turn-over' in j(Av). This slight 'turn- 
over' of 7, seen in Fig. 6, around Aw ~ 10 km s could 
be an indication of a significant contribution from warm 
gas with temperature < 10 3 K. 

In contrast, we do not see a similar converging trend 
for small values of Av end of this plot. Lazarian & 
Pogosyan (2000) suggests that typically for cold Hi Av ~ 
2.6km s _1 would be enough to approach the asymp- 



totic value for 'thin' slices (our thinest channel is 1.65 
km s _1 wide). Together with the observed turn-over 
at Av ~ 10, the lack of convergence for small Av sug- 
gests that Hi in the MB is primarily in the warm neutral 
phase. However, absorption observations of a number 
of continuum sources throughout the MB made by Kob- 
ulnicky & Dickey (1999), have revealed the presence of 
cold (< 20K) clouds. This cold phase is probably only 
a small fraction of the MB ISM and as such, the mech- 
anisms responsible for preserving the cold component in 
equilibrium with the warmer ambient ISM remain to be 
explained. 

4.4. Spatial power as a function of azimuthal angle 

By examining the power as a function of azimuthal 
angle, we can probe for any structure anisotropy of the 
MB Hi dataset. Any such anisotropies will manifest in 
the Fourier plane as significant structure as a function 
of phase, with excess power aligned aproximately in the 
direction perpendicular to that of the structure bias in 
the un-transformed image. Such azimuthal anisotropy 
could result from directionally-specific peturbations, for 
example, the effects of a magnetic field (Esquivel et al. 
2003), which may act to organise the ISM by shaping the 
ionised or dusty component. 

Fig.[S]shows the mean power as a function of azimuthal 
angle across spatial scales 100-600A of the Fourier trans- 
form for the Hi integrated intensity distribution of each 
of the four subregions. We find no evidence for any az- 
imuthal anisotropy in the Fourier plane, despite the fil- 
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Fig. 6. — Variation of the index 7 (P(k)ock 7 ) as a function of velocity for all four regions (marked in upper left) is plotted here as the 
solid line with error bars (left axis). The total integrated Hi brightness is over plotted as the dotted line (right axis). The shaded regions 
are those over which the power index is relatively constant and are the subjects for the study of the variation of spatial power-law index 

with integrated velocity thickness in Section 14.31 The velocity intervals shown here are those for which 7~0 i.e. those channels which are 
not dominated by noise. The horizontal dot-dashed line marks the RMS Tb for each subregion, as measured in line-free parts of the cube. 



amentary and apparently east-west orientation appear- 
ance of the Hi distribution in the southern regions shown 
in Fig. n 

It is likely that the dynamic range of any structure in 
the Hi images is insufficient to manifest as features in 
the Fourier transform: for such a feature to be distinct 
it is necessary to have a power at least 2-3 orders of 
magnitude larger than the surrounding structure. This 
is equivalent to an integrated temperature difference of 
10-30 Kkm s _1 . The apparent East-West structure of 
Fig. H shows an integrated brightness range significantly 
smaller than this, typically only ~2-5Kkm s^ 1 less than 
the ambient and less structured Hi. In addition, these 
structures are not particularly well defined, which will 
dilute the overall effect. 

5. THE SPECTRAL CORRELATION FUNCTION ANALYSES 

The motivation for the use of the SCF here is intended 
as a benchmark, by which we can compare the outcomes 
of the SCF and the SPS on the same dataset, and to test 
the response of the SCF on a purely tidally-generated 
system and in comparision to previous SCF studies of 
the LMC (Padoan et al, 2001). 

5.1. Derivation of the SCF 

The algorithm used here to calculate the SCF is that 
employed by Padoan et al. (2001) for the study of the Hi 
intensity distribution in the LMC. The SCF is defined 



in the following way (see also Rosolowsky et al. 1999 
Padoan, Rosolowsky, & Goodman, 2001): 



S (r,Ar) , 
with the numerator being defined by: 



5 (r,Ar) = 1- 



' j: v {T(r,v)-T(r + Ar,v)}2 
E v T(r,v)2 + j: v T(r + Ar,vy 



The noise term in the denominator is defined for the 
subject spectrum, T(r,v), by: 

5 ,Ar(r) = 1 - 



Q(r) 

where Q(r) is the spectrum 'quality' and is given by: 



(2) 



and R is the RMS in signal free parts of the cube. 

We aim to generate maps of the spatial variation of 
5*o(r,Ar), for a number of regions in the Hi MB dataset. 
The algorithm used is outlined below: 

1. Dataset Selection. A square region of a side with N 
pixels is defined within the data cube. The SCF is 
calculated iteratively using every spectrum in this 
region in turn as the subject spectrum. 
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Fig. 7. — Modification to the SPS index as a result of successively thicker velocity averaging for three of the four subrcgions shown in 
Fig. The SPS index of the North East region is quite distinct from the two southern regions and changes more rapidly at narrower 
velocity intervals than does the southern regions. 



2. Calculate Sq(t, Ar) map. A new larger region of 
an area of (2N — l) 2 is denned with the subject 
spectrum at the centre. The S*o(r,Ar) for each 
spectrum in this region is calculated relative to the 
centre subject pixel. This process forms a two di- 
mensional map, where the value at each point is 
proportional to the similarity of the spectrum to 
that of the centre spectrum. The value for the cen- 
tre pixel in this SCF map is unity and all other 
pixel values are <1. This map is stored for subse- 
quent averaging. 

3. Iteration and averaging. Step 2 from above is re- 
peated for each pixel in the initial N x N region 
of interest. Ultimately, a total of N 2 maps are ob- 
tained, each having an area of (2N — l) 2 . The final 
result is determined by finding the mean for all of 
the calculated SCF maps: S(Ar)= (S , (r,Ar)) r 

This approach requires a 'buffer' of spectra, at least 
N pixels wide around the region of interest. Without 
this buffer, there is not a uniform sampling population 
for each pixel of the SCF map and edge effects become 
important. The numbering in Fig. [3] shows the grid ref- 
erence of the centre square region defining the subject 
spectra (step 1 above). 5o(r,Ar) is calculated for all 
spectra within each box relative to all spectra within 20 
pixels (step 2 above) and then averaged over all 21 2 pix- 
els in that map (step 3). Overall, we obtain 8 different 
maps of So(Ar) from this dataset. 



5.2. Results of the SCF analyses 

The results for the SCF analysis are shown Fig. |5J The 
left-hand panels show the eight SCF maps, correspond- 
ing to each of the labeled regions in Fig. [3J Note that 
the maps contain significant structure, so that simply 
calculating So as a function of |Ar| = \J x 1 + y 2 (i.e. 
the azimuthally-averaged approach taken by Padoan et 
al. 2001) does not characterise the Sb(Ar) appropri- 
ately. As such, the right-hand panels of Fig. [3] show the 
behaviour of So(|Ar|) along the major and minor axes 
of the elliptical contours. The map is clearly anisotropic 
and indicates some structure bias present in the organ- 
isation of the gas in the MB. In particular, the major 
axes of the ellipses appear to be generally aligned with 
direction of filamentary structure in the MB. 

In order to interpret the SCF maps, we have fitted 
an elliptical power-law function to the distribution. We 
modeled the typical So surface using a function with pa- 
rameters Cf. 

f(x,y) = c {x 2 +y 2 )W 2 , where (3) 

f3 = c\ cos 2 ((/) — cs) + C2 sin 2 (</> — C3) , Ci > C2, and 

4> = arctan(y/x) 

This function characterises the So surface with two sep- 
arate power-laws falling off in orthogonal directions from 
Ar = with indices c\ (major axis) and C2 (minor axis) 
respectively. The angle C3 represents the orientation of 
these axes with respect to the coordinate axes. The data 
So(x,y) are weighted inversely proportional to their dis- 
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Fig. 8. — The variation of power in the Fourier transform of an integrated intensity map for each of the four regions is plotted here as 
diamonds. The log of Power is shown as a function of azimuthal angle. The dotted circle shows the mean power for this baseline interval. 
Error bars of one standard deviation show that there is no convincing consistent deviation from the dotted circle, and no indication of 
coherent structure in the Fourier transformed data for any of the four regions. 



tance from £0(0,0). The singularity at /(0, 0) for nega- 
tive values of c\ and C2 is ignored. We fit the function / 
to So for every position that is completely sampled. 

Padoan, Goodman, & Juvela (2003) point out a rela- 
tionship between the power-law scaling of So(| Ar|) and 
the line width of the spectra considered. We find no such 
correlation present in the MB Hi data, however, the 13 CO 
molecular cloud data used by Padoan, Goodman, & Ju- 
vela (2003) are systems with well defined size-line width 
relationships and where such line profiles can typically 
be described with only a single component. In contrast, 
the Hi data used here are complex, and show multiple 
velocity components for any line-of-sight (e.g. McGee & 
Newton, 1986; Muller et al. 2003) 

There is evidence for a correlation between the power- 
law index of the SCF and the integrated intensity. This is 
most notable when step 3 of the SCF processing is omit- 



ted and the major and minor axes are derived by fitting 
without averaging. Fig. ^] shows this correlation with 
respect to the major axis index, c\ , for all points outside 
the 21 pixel buffer at the edge of the map. The correla- 
tion is not outstanding, but appears to obey a power-law 
such that |ci| oc (f TAdv) 5±01 . Padoan, Rosolowsky, 
& Goodman (2001) note that higher density gas appears 
to be traced by a steeper power-law slope. However, this 
effect may be due to different spatial distributions for 
different molecular tracers. 

6. DISCUSSION OF STATISTICAL RESULTS 

6.1. Spatial Power Spectrum 

We find that the SPS of the Hi intensity distribution in 
the MB can be well described by a single power-law func- 
tion over a wide range of spatial scales (29 pc to at least 
2 kpc). As in the case of the SMC (Stanimirovic et. al, 
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Fig. 9.— Results of SCF analysis on the Hi in the MB for grids (2,2)-(5,2) (left panel; see Fig.0 and for grids (6,2)-(9,2) (right panel). 
The left part of each panel shows the map of So(Ar) according to the SCF algorithm in the text. The right part of each panel shows radial 
plots of S"o(| Ar|). The two cuts indicated in the left panel (solid and dashed line) appear in the right panel for the major (black) and minor 
(gray) axes of the elliptical So map respectively. We see that without exception, the ellipses are aligned with the major axis (i.e. the more 
slowly changing axis) approximately East- West. Each panel is annotated to indicate the corresponding Hi subset defined on Fig. [3] 



1999a, b), no significant changes to the power-law shape 
of the SPS were found when varying heliocentric velocity 
or the thickness of velocity channels, though the power- 
law index varies significantly. This indicates that the 
MB is likely to have a large line-of-sight depth. In con- 
trast, the SPS derived from the Hi intensity distribution 
in the LMC showed a break at the spatial scale equiv- 
alent to ^100 pc (Elmegreen, Kim & Staveley-Smith, 
2001) which was interpreted as evidence for the finite 
thickness of the LMC disk. However, the SMC and the 
MB have large-scale morphologies which are distinctly 
different from the LMC, namely that they are not ap- 
proximately face-on systems and have a substantial line- 
of-sight depth. 

There is a clear variation of ~/(Av) over the range of 
velocity widths considered, and the NE section of the 
analysis has qualitatively different variation than the SW 
and SE sections. Lazarian & Pogosyan (2000) relate the 
measured variations of 7(Aw) to physical properties of 
the turbulence, a particular strength of SPS analysis. In 
particular, j(Av) can be related to the indices of the 
three-dimensional power spectra of the density and ve- 
locity fluctuations. If the density fluctuation power spec- 
trum is a power law with index of n and the velocity 



fluctuation power spectrum has an index of fi the rela- 
tionships between 7, n and fi are given in Table ^ for 
n > —3 (Shallow 3-D density) and n < —3 (Steep 3-D 
density) . 

Following the description of Lazarian & Pogosyan 
(2000): if we assume that the smallest Aw traces 'thin' 
velocity slices and that the largest Av corresponds to 
'thick' velocity slices, we can estimate the factor by which 
the velocity fluctuations modify the SPS density spec- 
trum. The derived values for the slopes of the density 
and velocity spectra are listed in Table [21 These values 
are given as upper and lower limits since asymptotic so- 
lutions were not fully met, especially at the lower Av 
end. 

Similarly, if we assume that the thickest velocity 
slices used in the study of Section 14.31 are probing the 
pure density-dominated regime of Lazarian & Pogosyan 
(2000), then the 3-D density index for the SW region has 
a lower limit of n rj —3.4, indicating that this region is 
in the velocity-dominated regime. Using the SPS index 
for the thinnest velocity slice of this region leads to a 
lower limit estimate of the velocity fluctuations /1 ~ —4 
(Table 2). Similarly for the NE region, the SPS index 
for thickest velocity slice has a lower limit of n » —2.9, 
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Shallow 3-D density Steep 3-D density 


Thin slice (Av small) 
Thick slice (Av large) 
Very Thick Slice (Ad very large) 


n - n/2 - 3/2 -9/2 - fi/2 
n -3/2 + fi/2 
n n 



TABLE 1 

The value of the observed power-law index 7 FOR. different values of Ad in different turbulent regimes, n IS THE STATIC 

SPATIAL POWER LAW INDEX OF THE TURBULENCE WHICH IS MODIFIED BY VELOCITY FLUCTUATIONS WITH A SLOPE OF y. (LAZARIAN & 

POGOSYAN, 2000). 



indicating that this region is in the density-dominated 
regime and that velocity fluctuations are small in this 
region. The derived 3-D velocity index has an upper 
limit of yu ~ —4.5 and is extremely steep, indicating an 
overwhelming deficiency of high- velocity fluctuations. 

The steepening of the SPS slope with increasing in- 
tegrated velocity thickness, as predicted by Lazarian & 
Pogosyan (2000) , was also measured from Hi in the SMC 
(Stanimirovic & Lazarian 2001) and for a test region in 
the Galaxy (Dickey et al. 2001). Table 2 shows that the 
slopes for the southern two regions in the MB are gen- 
erally similar to those derived for the SMC which has a 
n = —3.3 and \i = —3.4 (Stanimirovic & Lazarian 2001). 
Both the SMC and the MB appear to have power-law 
indices that are more shallow that what is predicted for 
a Kolmogorov turbulent spectrum. For the case of Hi 
in the Galaxy a shift of the SPS slope from about —3 
to —4 after velocity integration was noticed by Dickey 
et al. (2001). Other studies of the LMC (Elmegreen, 
Kim & Staveley-Smith 2001) and the Magellanic Stream 
(Stanimirovic et al. 2002) also show apparent steepening 
of the power spectra with increasing integrated velocity 
thickness. 



A different interpretation of the power-law SPS in 
the case of the SMC was provided by Goldman (2000) 
whereby a featureless SPS was interpreted as a signature 
of turbulence being driven on spatial scales larger than 
the size of the SMC. This turbulence was induced by lo- 
cal instabilities due to tidal interactions in the Magellanic 
System. This approach also defines the relationship be- 
tween the 2-d intensity and 3-d density statistics, where 
the Power-law index of intensity fluctuations, 7, is re- 
lated to the power-law index of density fluctuations, n, 
by 7 = n + 1 . In the MB case studied here, this approach 
will result in n <~ —4 for the Southern two regions, and 
n ~ —3.2 for the North-East region. 

The second important result from these studies of the 
MB, is that the slopes of the SPS derived for the northern 
and southern regions appear to be strikingly different: 
the northern part shows a distinctly shallower (although 
still monotonic) power spectrum, indicating a relative 
lack of large-scale power (or equivalently, an excess of 
small scale power). Furthermore, these results could be 
interpreted as evidence for a warm (~10 3 K) component. 
This is slightly perplexing, given the results of absorption 
detections which indicate the presence for a very cool 
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iThin ^Thick (= n) /i 3-D Velocity 

(lower limit) (upper limit) (upper limit) 
3-D Density 

South East -3.05±0.06 -3.4±0.2 -3.6±0.3 

South West -2.94±0.04 -3.4±0.2 -3.9±0.2 

North East -2.24±0.03 -2.9±0.2 -4.5±0.4 



TABLE 2 

SPS INDICES FOR THICKEST AND THINNEST VELOCITY SLICES FOR THE THREE SAMPLED REGIONS IN THE MB (SE, SW AND NE REGIONS OF 
FlG.0. THE DEGREE BY WHICH TWO-DIMENSIONAL VELOCITY FLUCTUATIONS AFFECT THE SPS INDEX ARE SHOWN, ALONG WITH THE 
ESTIMATED 3-D VELOCITY INDEX (SEE LAZARIAN & POGOSYAN, 2000 AND TABLE . 



(~20 K) component. 

Muller et al. (2003) have observed that at least two 
distinctly different Hi velocity components exist through- 
out the MB, being centered on about 150 and 190 km s _1 
(see also Fig. Furthermore, the two components are 
distinct in declination: the higher velocity component is 
predominantly at northerly declinations and is sampled 
by the NE and NW regions, while the lower velocity com- 
ponent is sampled by the SE and SW regions. 

Numerical simulations (Gardiner, Sawa & Noguchi, 
1994) suggest that the two velocity components could be 
associated with arms extending from the SMC. In this 
model, the more distant arm is projected as a higher ve- 
locity and slightly northern feature, while the nearer arm 
correlates with lower velocity component. The nearer 
arm is clearly contiguous with the SMC and LMC may 
be more properly regarded as 'The Magellanic Bridge'. 
The distances to these two arms have been predicted by 
Gardiner, Sawa & Noguchi (1994) to be centered at 50 
and 70 kpc, with an expected line-of-sight depth for both 
entities of ~5-10 kpc. We propose that as an explanation 
for the differences in the behavior of the power indices of 
the northern and southern regions in Figs.[S]&[7| these 
regions encompass two separate and distinct arms of the 
SMC. We view the arms projected approximately on top 
of one another, whereas these features are not well spa- 
tially connected at all. 

6.2. Spectral Correlation Function 

Like the SPS, the SCF shows a logarithmic dependence 
of scale which can be parameterised with a single compo- 
nent. We do not expect the slopes of the SCF and SPS 
to be consistent, since the SPS probes structure in excess 
of the mean, large scale power (i.e. that represented by 
the zero lag), whereas the SCF algorithm, as used here 
and by Padoan et al. (2000), does not factor for a mean 
brightness level which will otherwise tend to flatten the 
SCF slope. Interestingly, the SCF slopes derived here 
for the MB dataset are largely consistent with those de- 
rived from the Hi in the LMC (S (Ar) 0.15 - -0.4; 

Padoan et al, 2001). 

The lack of radial symmetry in the maps of So(Ar) 
is itself an interesting result. It appears then that the 
SCF is more successful in detecting structure trends over 
a smaller dynamic range than the approach using the 
Fourier transform and the SPS in Section 

The obvious structure in the approximate East- West 
direction (i.e. the line joining the SMC and LMC) indi- 
cates a more persistent similarity of spectra along that 



direction. What we see here is that the rate of change of 
the entire spectra (rather than a single velocity channel, 
or integrated brightness) across the entire sampled ve- 
locity range is more gradual in the East- West direction 
than in the North-South direction. 

We interpret the East- West 5o(Ar) structure as a re- 
sult of the tidal stretching imposed by the LMC-SMC in- 
teraction, although in this case the morphology of the Hi 
is already visually suggestive of this process. We should 
also bear in mind that the brightest parts of the dataset 
can dominate the SCF results, biasing steeper power-law 
indices and the overall structure of the SCF maps. A 
floating temperature scaling factor in the SCF algorithm 
would probably go some way to removing this bias. 

We note that there is no evidence for a break in the 
Bridge SCF power law, as was found for the LMC at 
^100 pc, although these two objects are widely sepa- 
rated in space and there is no outstanding reason that 
they should share a common morphology at that scale. 
However, since a typical width of the structures in the 
filaments is ~ 100 pc in the plane of the sky, and we 
do not see a departure from a power law on this scale, 
this may be another indication of significant line-of-sight 
depth in the Bridge. We see in Fig. that there may 
be some evidence for a consistent break in the power-law 
for both the major and the minor axes, at ~300 — 500 
pc. However, this departure occurs at scales which are 
comparable to the maximum size of the tested areas (i.e. 
^500 pc). The degree to which the boundary and size 
of the sample affects the SCF has not been well studied, 
and it is probable that the departure is influenced by the 
edge of the sample region. As such, any interpretation 
for the reasons of the observed departure at this scales 
should be made with caution. It is not appropriate to 
simply expand the spatial size of the sampled areas since 
this would sample data which is excessively dominated 
by noise and the SCF algorithm will be affected unpre- 
dictably. 

In this study, we did not attempt to study the be- 
haviour of the SCF in a similar way to that used by the 
velocity-component analysis (VCA), namely, measuring 
the modification to the SCF as a function of varying ve- 
locity integrated width. As the SCF operates in the im- 
age domain, we expect that by averaging across the ve- 
locity fluctuating component the SCF will flatten some- 
what, and that there will be some correlation of the SCF 
slope with velocity integrated width. However, to con- 
duct this test appropriately, a study with a well charac- 
terized dataset should be used, rather than the complex 
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Hi profiles of Magellanic Bridge dataset. 

7. SUMMARY 

We have attempted two statistical analyses of the fila- 
mentary Hi in the Magellanic Bridge. Together with the 
fact that there also appears to be an identifiable disconti- 
nuity in the velocity distribution of the Hi in the MB, the 
most significant result from the calculation of the spatial 
power spectrum (SPS) here has provided support for a 
scenario suggested by numerical simulations, where the 
southern and northern parts of the Bridge represent the 
projection of two distinct arms of gas emanating from the 
Small Magellanic Cloud (SMC). These findings suggest 
that we now need to re-assess the current interpretation 
of the Bridge as a single filamentary feature. 

Similarly to the SMC, the power spectra throughout 
the MB are well described as a power law with a sin- 
gle power component. In particular, the southern region 
which contains brighter and more turbulent Hi, have 3- 
d density and velocity indices similar to that which was 
found for the SMC. The more tenuous and higher veloc- 
ity (i.e. northern) parts show a significantly shallower 
SPS, with the 3-d density slope of —2.9. This is indica- 
tive of relatively quiescent gas. The 3-d velocity index 
(—4.5) of this region is significantly steeper that for the 
case of a Kolmogorov turbulent flow. 

We find that the a velocity component analysis (VCA) 
of the Hi in the MB behaves, on the whole, according to 
predictions by Lazarian & Pogosyan (2000), where mod- 
ifications to the power spectra by velocity fluctuations 
can be removed by integrating over a velocity window. 
This has again highlighted the dissimilarities between the 
northern and southern parts of the MB. The contribu- 
tions by velocity fluctuations to the southern part of the 
MB are again similar to those in the SMC, whereas the 
SPS from the northern part of the MB shows evidence 
for a very steep three-dimensional velocity index and a 
lack of rapid velocity fluctuations. 

A study of the azimuthal isotropy of the Fourier trans- 
form of the Hi dataset does not show any such East-West 
structure which is apparent in the Hi peak and integrated 
intensity maps. It is likely that the dynamic range of the 
Fourier transformed dataset is insufficient to distinguish 



any such structure tends from the ambient Hi. Con- 
versely, the spectral correlation function (SCF), due to 
the nature of the algorithm, is much more successful in 
indicating more subtle, local trends in structure. 

The analysis using the SCF has shown more quantita- 
tively the effects of tidal stretching on this dataset. It 
confirms that the observed Hi spectra in this region vary 
more slowly along the approximate direction of the tidal 
stretching. Although limitations of the application of 
the SCF has confined the region sampled to the bright- 
est parts of the MB, the SCF shows enormous potential 
as a tool for more easily parameterising suspect tidal 
features, insofar as the direction and extent of the tidal 
perturbation. 

Although the behavior of the SCF cannot be directly 
related to the physical properties of the underlying tur- 
bulence, the SCF provides some provocative results. We 
see that, like for the SPS, the SCF shows a power-law 
dependence with indices compatible with previous stud- 
ies of the LMC Hi datset. However, the results from the 
SCF analysis also hint at a characteristic scale length, 
which is not observed in the SPS of the same region and 
are not at scales consistent with the studies of the LMC. 
It is important to bear in mind that the departure ob- 
served here is at the limits of the sampled scale range 
and any interpretation should be made carefully. How- 
ever, we can conclude that the SCF and the SPS do not 
appear to measure structure variations in the same way. 
Specifically, the SCF appears to be much more sensitive 
to low-power and small scale variations in structure. It 
will be necessary to explicitly test and compare behav- 
iors of the SCF and SPS to different turbulence types 
before the full application of the SCF can be recognized. 
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